www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第四讲 随机数的生成及随机变量抽样/liti7.m

    mm=10000
Randnum=unifrnd(0,4,1,mm);
xRandnum=zeros(1,mm);yRandnum=zeros(1,mm);
for ii=1:mm
    if Randnum(1,ii)<=1
       xRandnum(1,ii)=0;
       yRandnum(1,ii)=Randnum(1,ii);
       else if  Randnum(1,ii)<=2
            xRandnum(1,ii)=Randnum(1,ii)-1;
            yRandnum(1,ii)=1;
                else if  Randnum(1,ii)<=3
                         xRandnum(1,ii)=1;
                          yRandnum(1,ii)=1-(Randnum(1,ii)-2);
                      else     
                          xRandnum(1,ii)=1-(Randnum(1,ii)-3);
                          yRandnum(1,ii)=0;
                      end 
            end
      end     
 end      
xRandnum
yRandnum
 [Y,JJ] =sort(xRandnum)
plot(xRandnum(JJ),yRandnum(JJ),'.')